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ABSTRACT 

Linear fluid equations are used to estimate wave train properties 

of strong collisionless shocks. Fast shocks exhibit several dispersion 

changes with increasing Mach number. For perpendicular propagation 

into a finite-0 plasma, an ion cyclotron radius trailing wave train 

2 

exists only for M„ <2. Oblique fast shocks have a leading ion 
inertia wave train if M. < /M , /M cos 6/2 and a trailing electron 

A T - 

inertia wave train if M > /mH/M cos 0/2. If the downstream sound 
speed exceeds the flow speed, linear wave theory predicts a trailing 
ion acoustic structure which probably resides within the magnetic 
shock. For a turbulent shock model in which an effective electron-ion 
"collision" frequency exceeds the lower hybrid frequency, ions decouple 
from the magnetic field; the shock wave train now trails with electron 
inertia and electron gyroradius lengths. Comparisons of this turbulent 
model and observations on the earth's bow shock are made. 
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1.0 Introduction 


In the fluid theory of collisionless shocks the transition of 

macroscopic plasma quantities from a steady upstream to downstream flow 

«• 

state is described by a spatially oscillatory wave train. Oscillation 
scale lengths are characterized by the natural plasma dispersion lengths, 
familiar from the theoi'y of linear wave propagation [Stringer, 1963]. 
Laminar solutions have been applied primarily to weak collisionless 
shocks, partly because of analytical simplicity and partly because of 
fundamental difficulties inherent in the fluid formulation [Sagdeev, 

1966; Cavaliere and Englemann, 1967]. This paper considers extension 
of fluid theory into the strong shock regime concentrating exclusively 
on changes in the wave train dispersive structure. Such restriction 
permits utilization of linear concepts to obtain qualitative and semi- 
quantitative estimates of the shock structure. 

Fluid shock theories have evolved on two distinct and clearly 
separable levels: dissipative shocks in which Coulomb interactions 

control the structure, and dispersive shocks in which collisionless 
plasma properties dominate. For the present purpose the dissipative 
limit is instructive since here the structure of strong shocks is 
reasonably well understood. The transition between allowed Rankine- 
Ilugoniot states occurs over dissipative scale lengths such as the 
resistive, viscous, or thermal diffusion lengths [Anderson, 1963]. Since 
dissipation rates vary inversely with thickness, strong shocks, which 
require intense overall randominzation , often possess a multiply 
dissipative structure. For example, for fast perpendicular shocks in 



which the downstream sound speed exceeds the flow speed, a condition 

required by the Rankine-Hugoniot relations for upstream Mach numbers 

exceeding about two, the fluid-magnetic, field resistive coupling fails 

to provide the necessary dissipation. A thin viscous discontinuity or 

subshock in which the fluid and magnetic field are decoupled and intense 

* 

viscous dissipation is generated develops within the broader resistive 
structure [Marshall, 1955; see Fig. 1]. Similar dissipation discontinuities 
occur in a variety of other fast and slow shock flows [Coroniti, 1970]. 

In the collisionless limit dissipation processes are assumed 
sufficiently weak to permit wave dispersion domination of the shock 
structure. The consistent approximation of weak shocks reduces the 
complications from nonlinear fluid behavior and allows some confidence 
in the validity of wave train concepts based on the linear wave theory. 

Weal; collisionless dissipation, such as results from micro-instabilities 
and three wave decays [Sagdeev, 1966] is usually considered to be 
driven by the wave train spatial gradients ; the total dissipation is 
accumulated over many wave train oscillations. Hence weak dispersive 
shock models have generally avoided discontinuities within the wave 
train structure. 

Theoretical considerations of strong collisionless shocks have 
previously stressed the development of fully turbulent flows in which 
fluid effects are unimportant [Fishman et al., 1960; Camac et al., 1962; 
Kennel 'and Sagdeev, 1967; Tidman, 1967], Although dissipation again 
primarily controls the structure, the nonlinear turbulent interactions 
are so complicated that often no assurance of a complete shock transition 
in the Rankine-Hugoniot sense is possible; i.e., the downstream state 
remains unspecified. Therefore, given these difficulties with fully 



turbulent shocks, the possibility exists that fluid theory might prove 

a useful guide to at least some aspects of strong shock structure. 

Alternatively, strong turbulent dissipation has the tendency of making 

the plasma behave like a fluid on macroscopic scale lengths; hence some 

remnants of fluid structure might be expected for strong shocks. 

* 

In analogy with multiply dissipative collisional shocks, strong 
collisionless fluid shocks possibly contain dispersion discontinuities 
or multiply dispersive wave trains. Strong shocks, however, require 
strong dissipation. Since plasma turbulent dissipation should be tightly 
coupled to wave train dispersion gradients, the separation between 
dispersion and dissipation, which permits great analytical simplification 
for weak shocks, is considerably more suspect for strong shocks. Recent 
observations on the earth's collisionless bow. shock, however, provide 
some encouragement for this approach. Fredricks et al. [1970] report 
that the magnetic shock is often characterized by an electron inertia 
length wave train- like structure; furthermore, strong electric field 
wave turbulence, thought to provide considerable ion heating, is coupled 
to the magnetic field gradients. Hence fluid dispersion remains even 
in the presence of rapid turbulent dissipation. 

The analysis of the nonlinear wave train structure including the 
self-consistent turbulent dissipation is not yet feasible because of 
both computational complexity and the absence of a complete theoretical 
formulation for strong plasma turbulence (see Kellogg [1964] and 
Bardotti et al. [1967] for some initial efforts). Therefore, the 
philosophy of thispaper is to continue utilization of linear fluid 
arguments to obtain qualitative estimates on the nonlinear behavior of 
strong collisionless shocks. In the vicinity of the shock Rankine-Hugoniot 



states, wave train amplitudes are small; here the two fluid theory of 
linear waves can be used to predict dispersion scale lengths and the 
sign of dispersion, i.e., whether the wave train leads or trails. A 
concomitant approach is to derive a linearized differential equation for 
the wave train whose coefficients are then evaluated according to the 
Rankine-Hugoniot relations. This method is somewhat less useful, however, 
since a consistent set of fluid equations in differential form which are 
valid over a wide range of dispersion scale lengths is difficult to 
formulate. Although the nonlinear structure of the shock center is not 
realized by these techniques, numerical computations, for restricted 
shock parameters, have substantiated the general features of the w ave 
train structure predicted by linear theory [Cavaliere and Englemann, 1967] 
Section 2.0 reviews the dispersion limitation to the steepening of 
a finite amplitude compression pulse. Changes in the perpendicular 
fast shock wave train arising from multiply dispersive wave propagation 
are considered in terms of the ion acoustic subshock proposed by the 
Texas group (Robson, 1969) and others (Paul, 1969), 

Sections 3.0 and 4.0 discuss fast shock wave trains from several 
points of view. In Section 3.0 an ion gyroradius length wave train is 
derived for a finite-3 perpendicular shock. For Mach numbers less than 
two, the predicted shock thickness becomes less than the ion gyroradius, 
thus violating the small gyroradius expansion of the fluid equations; 
a speculation on the structure of stronger shocks is presented. The 
restriction of perpendicular propagation is relaxed in Section 4.0. 

Oblique fast shocks exhibit a dispersion change from a leading ion 
inertia to a trailing electron inertia wave train for upstream Alfven 
Mach numbers exceeding /M /M cos0/2. Finally the ion acoustic subshock 
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is reconsidered for oblique propagation; restrictions on its occurrence 
in terras of propagation angle, plasma 3, and kinetic theory effects 
are discussed. 

Section 5.0 attempts to simulate the possible effect of a strong 
turbulent resistivity on the wave train dispersion. If the anomalous 
collision frequency exceeds the lower hybrid frequency, ions are decoupled 
from the magnetic field, thus removing ion inertia dispersion even for 
oblique propagation. The strong shock magnetic structure is then 
controlled by electron inertia and electron gyroradius dispersion, both 
of which produce trailing wave trains. 

Section 6.0 considers the observations on the earth's bow shock in 
terms of the previous discussion, and concludes that turbulent resistive 
decoupling of ions is probably required to explain the electron inertia 
magnetic wave train. A speculation on the structure of slow shocks 
predicted to occur in neutral sheet flows is presented. 



- 8 - 


2.0 Steepening Limitation of- Coll isionless Shocks 

Finite amplitude, long wavelength compression pulses steepen to form 

shocks [Kantrowitz and Petschek, I960]. The shock structure is determined 

by whatever processes limit fluid gradients from becoming arbitrarily 

large. In collisional shocks, steepening ceases when the pulse width 

becomes comparable to dissipation scale lengths. In collisionless plasmas 

steepening limitation arises from linear wave dispersion, which is reviewed 

here following the arguments of Sagdeev [1966], 

As the pulse propagates, nonlinear terms in the fluid equations 

generate harmonics of the fundamental pulse wave frequency (o.O- or wave 

number (k) . If the harmonic-waves also satisfy the same non-dispersive 

dispersion relation, phase speed to/k = constant, as the pulse fundamental 

mode, the pulse propagates as an entity and steepening continues. Hence 

whenever the small amplitude dispersion relation predicts non-dispersive 

propagation over some range of w and k [as in Fig. 2 for w/k = 

2 2 1/2 

CC A + Cg ) and w/k = Cg ] , compression pulses steepen. As steepening 

proceeds, however, harmonic modes with’ wavelengths comparable to plasma 

dispersion lengths are excited. Typical dispersion lengths are the ion 

and electron inertial lengths, C/a>p , the respective particles' gyro- 

radii, R + = C + /|n + [ , and the particles Debye lengths, = C + /u)p 

“ ± ~ ± 

The i signs refer to ions and electrons, respectively; = ± eB/M + C 

is the gyro- frequency , B is the magnetic field strength, C is the 

2 1 /? 

velocity of light; w p = (4ttNc /M + ) is the specie plasma frequency, 

‘ + - 

N is the number density, e the electronic charge, M the specie mass; 

+ 1/2 ± 

C + = (T /M + ) is the particle thermal speed, and T is the 

temperature in energy units. Gaussian units are used throughout. 
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To illustrate the dispersion limitation of steepening, consider th 

particular example of the linear fast wave propagating perpendicular to 

. ± ± 2 
the. magnetic field. The dispersion relation when 3 = 8 ttNT /B < 1 , 

T /T < 1 , and Wp > |f2 + Q is approximately [Stringer, 1963] 

-Vt- ’ 

T (1 + k \ ) 


w A 

2 2 2 2 
k 1 + (k/cVtOp ) 


S 2 

+ C 

' , 2,2 L + 

1 + k X D 


( 2 . 1 ) 


2 1/2 

and is sketched in Figs. 2a and 2b i = (B /4 ttNM + ) ' is the Alfven 

1/2 

speed and C = (T /M ) is the ion acoustic speed. Note that the 

O * “ 

fast wave is multiply dispersive with both C/( 0 p and dispersion 

lengths. From (2.1) the longest dispersion length attained by the 

harmonic waves is C/cJp ; electron inertia decouples the fast wave from 

. 6 , 

the magnetic field and decreases the phase speed below the fast hydro- 

2 2 1/2 

magnetic speed, w/k = (C^ + C<, ) = C p . 

Harmonic waves with kC/o^, n, 1 , therefore , propagate slower than the 
pulse and trail behind. Steepening limitation and the formation of a 
steady state flow is now possible if the compression energy of steepening 
is removed by wave energy convected out of the pulse. If thermal corrections 
to (2.1) are neglected, (T .<< 1 , the steady pulse is the magnetic soliton 
of thickness C/w p [Ad lam and Allen, 1958; Sagdeev, 1958; Gardner et al., 
1958; Davis et al., 1958]. If the plasma is weakly dissipative, this soliton 
is converted into a shock wave train [Sagdeev, 1966] consisting of trailing 
dispersively propagating waves which, in the shock frame, phase stand in 
the downstream flow. Since, by the steepening limitation arguments, the 
wave train must convect energy away from the shock, the group velocity of 
the trailing waves must be less than the phase velocity, 3w/3k < to/k . 
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From (2.1) note that % pr 2^'2 — 2 — < F ^ 3 « 1 . 

(1 + k C /iOp_ ) 

Asymptotically downstream the wave train amplitude, having been 
weakly damped by dissipation, is sufficiently small that the waves 
approximately obey the linear dispersion relation. Hence an estimate of 
the wave train oscillation scale length, L , can be obtained by equating 
the phase speed from (2.1) to the downstream flow speed predicted by 

the Rankine-Hugoniot relations. (In Fig. 2a this criterion is satisfied 
at the intersection of with the linear dispersion curve, provided, 

of course, that ,3w/3k < to/k at intersection.) Thermal corrections to 
(2.1) for this estimate are small if or < 1 - Performing 

the calculation yields 



where the downstream Mach number M„ = U~/C,. . The spatial structure 

*2 1 P 2 

of the magnetic wave train is also sketched in Fig. 2a. 

When kC/Wp » 1 , the fast wave increases its phase velocity to 
the ion acoustic speed, C g , and again propagates non-dispersively , 
w/k * Cg until Debye wavelengths are reached; for kX^ > 1 the wave 
slows to the ion thermal speed, C + . Therefore in the locally non-dispersive 
region, w/k - C Q , . further shock steepening can occur; indeed when 
U 9 < C , the shock must develop short scale structure at the Debye 

length since the wave train phase standing condition can only be satisfied 
by a trailing ion acoustic wave. Note from (2.1) that 3w/3k < w/k for 
this wave. On this basis the Texas group [Robson, 1963] has suggested 
that a general feature of moderately strong collisionless fast shocks is 
the presence of an ion acoustic subshock. Since the shortest length over 



which the magnetic field can change is C/w p >> , and also since the 


ion acoustic wave is almost electrostatic [Formisano and Kennel, 1969], 
the magnetic field is virtually unaffected by the ion acoustic wave train. 
The subshock acts primarily to slow and heat the ions, and therefore can 
be considered as the dispersion analogue of the dissipative viscous 

t 

subshock. 

The upstream flow conditions for which an ion acoustic subshock must 
occur are determined by the perpendicular shock Rankine-Hugoniot relations, 
Skipping the trivial algebra, CL /U„ as a function of the upstream Mach 

^ j O ^ 

2 2 Z 2 

number M„ and 6. = C„ /C, C 0 = 2P/p . is given by 


Ip and B 1 = Cg /C A C g = 2P/p , is given 


U, 


m; + 2 
F 1 . 


27 M, 


4IC - 1 
F 1 


(Mp^+2) 2 (1+3 1 ) J 


(2.5) 


Figure 3a is a graph of C„ /U„ against M with B. as a parameter; 

2 *1 

Fig. 3b sketches the critical Mach number, M . , for which C„ . = IL 

cn l ^2 4 

against B. . Shocks with M > 2.5 have ion acoustic subshocks for 
1 F 1 

any upstream B . The Rankine-Hugoniot relations, Eq. (2.3), actually 

yield an upper limit to the criteria U_ < C . The wave train will 

b 2 

generally have a nonlinear overshoot in the first oscillation [Sagdeev, 
1966]; hence in the leading edge U < Cg will occur for lower Mach 
numbers than given by (2.3). A sketch of an ion acoustic subshock 
occurring in the leading edge of a C/w p magnetic w'ave train is included 
in Fig. 2a. 

The oscillation scale length for the ion acoustic wave train can be 
estimated from the intersection with the linear dispersion curve. 


Since for reasonably strong shocks B 2 >> Bj_ , the dispersion curve differs 
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from that appropriate to upstream flow conditions; Figure 2b is a sketch 
of (2.1) for 3 2 > 1 . Setting = w 2 /k 2 = C" /(l + k^A 2 ) , the scale 

length is 


L * U - Ms 2 

l7TZp72 

b 2 


where NL 


becomes 


U 2 /C S 2 and X D ^ C S 2 /U P + 


(2.4) 

Substituting from (2.3) L 



/3 M 

*'l 




(4M 2 +1) (M 2 +2) 2 - [27Mp /(1+B.)] 
r l 1 j 1 

3(M 2 -1) (M 2 +2) 2 - [27Mp /(1+0 )] 
F 1 F 1 F 1 1 


(2.5) 


In summary, this section has established the following points on 
the utilization of linear fluid theory to determine wave train structure: 

a) Hydromagnetic compression pulses steepen in regions of non- 
dispersive linear wave propagation. 

b) Steepening is limited at pulse thicknesses comparable to dispersion 
scale lengths. Wave trains are formed by dispersively propagating waves 
which phase stand in the downstream flow and have 3w/3k < w/k , or in 

the upstream flow with 3co/3k > w/k . • ■ 

c) Wave train scale lengths can be estimated by matching linear w r ave 
phase speeds to flow speeds at the asymptotic Rankine-Hugoniot states. 

Linear dispersion curves appropriate to upstream or downstream Rankine- 
Hugoniot conditions must be used for these estimates. 

c) Since the linear dispersion relation is multiply dispersive, 

perpendicular fast shocks possess a magnetic wave train with C/ojp lengths, 

and for M_ > M , an ion acoustic subshock at Debye lengths. Because 
F, crit ’ 
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of the wave train nonlinear overshoot, the Rankine-Iiugoniot relations 
provide an upper limit on . Since it interacts primarily with 

ions leaving the magnetic field unaffected, the ion acoustic subshock is 
the dispersion analogue of the viscous subshock. 
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3.0 Tne Per pendic ular Ion C yclotron Rad ius W ave Train 

In the previous section the discussion of the perpendicular fast 
shock was restricted to consideration of only electron inertia and Debye 
length dispersion. As this shock steepens in a i finite-$ + plasma, however, 
the longest dispersion length encountered is the ion cyclotron radius 
(ICR). In Section 3.1 an ICR wave train differential equation is derived 
for the linear perturbed plasma response about the asymptotic stationary 
flows. For moderate Mach numbers the shock steepens to lengths shorter 
than R + thus passing beyond the validity of the fluid analysis. Kinetic 
theory effects and steepening to shorter scale lengths are discussed. 

3. 1 Derivation of the ICR Wave Train Differe ntial E quation 

To investigate the perpendicular ICR fast shock wave train, the 

analysis employs the Chew-Goldberger-Low hydromagnetic equations with 

first order ICR corrections as derived by MacMahon [1965], (A similar 

approach has been used by Kinsinger and Auer, 1969; Goldberg, 1970.) 

The range of validity for these equations requires 3 + >> M /M , assuring 

R >> C/w n , T + /T » (M /M )2 so that R >> R and Q /w n « 1 or 
- + 

R + >> , thus restricting consideration to quasi-neutral plasmas. In 

+ 

addition, only shock scale lengths- L >> R + can be treated, consistent 
with the first order ICR expansion of the fluid equations. The wave 
train structure is investigated only about the upstream and downstream 
stationary Rankine-Hugoniot flows, thus permitting a linearized analysis. 
This method examines the stability at the stationary points of the 
nonlinear- wave train differential equation; shock transition requires 
instability upstream and stability downstream. 
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The coordinate system in the steady shock frame, used throughout, 
this paper, is specified by shock propagation in the X-direction, along 
which plasma quantities vary, with the magnetic field, B z , in the 
2-direction. The perturbed time independent fluid equations are 


n , B 6 B 

pUdU + 6 P l J ~ = 

xx 4 tt 


d <5U 
dx 


, d 6V 

plJ <$V + 6 = rj — y - 

y xy dx 


6 pO) pU) 

pU 2 6U + UC6P, fl) + 6P (1) ) + pf 1 ' 1 6U + U —4— + ~ — 

xx ■* 2 2 


Ou 


, s,( 1) x e.(l) s ~ 

+ (oqp + <5q v ') • X + 

~ '-A 


U, B 6B 

1 Z 1 Z d dU 

4tt ;iU dx 


(3.1) 

(5.2) 


(5.3) 


U 6B + 
z 

B 6U 
z 

= 0 


(3.4) 

pU) 





6(U -f- 

* q“ (I) • i) 

= 0 

(3.5) 

6P (1) = 

XX 


p(0) + 

d i$V 


A 

2Q 

+ 

y 

dx 

(3.6) 


pP>) 

+ 



6P<« - 

xy 

2iT + 

d 6U 
dx 


(3.7) 



2P^ + 

a «p (1 > 


A 

X = 

A 

NH ft + 

xy 

dx 

(3.8) 


Following the notation of MacMahon [1965], superscripts (0) and (1) 
denote zero and first order quantities, respectively, in the small ICR 
expansion. 1 (U) denotes perpendicular (parallel) to the magnetic 
field. 6 denotes a perturbed variable; the unperturbed quantities are 
to be evaluated about either the upstream, subscript (1) , or downstream. 
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subscript (2), flow conditions. U is the flow velocity in the X-direction; 

Vy is the fluid velocity in the transverse Y-direction. P , appropriately 

sub- and super-scripted, denotes the pressure. 6q'*“^ * x and 6q*?^ • x 

'■'A 

are the first order ICR heat flows along X , and are retained since they 
contribute to the same order as the pressure tensor ICR terms in the 

energy equation (3.3); zero order heat flow has been neglected. Effective 

4 

viscosity coefficients, q , £ , and p = j q + C are included to 

simulate weak collisional or turbulent dissipation. 

Neglecting products of ICR and viscosity terms, the above equations 
reduce to a wave train differential equation for 6U 


2 2 
tp _ r *■ 

d oU F 


3 n 2 d iU . yU 

4 * P| (0)* ‘ax- - p |0), /p 


- 61! 


0 


(3.9) 


R + = (fj^ + /P^ + ^)^ 2 is the ion cyclotron radius, and 
2 m 1/2 

Cp = [C + (2F^ J /p)] is the fast hydromagnetic speed. Solutions 
of (3.9) in the form exp(Xx) are 


X = 


•yu 


3 p (0) + 

2 A + 


+ 

2 ~ 


u 2 - c f 2 ^ 1/2 


+ L 4 \t 


3 D (0) + 


/P 


(3.10) 


where terms proportional to p have been neglected. About the upstream 
flow U exceeds Cp by the shock evolutionary conditions [Kantrowitz 
and Petschek, 1966]; (3,10) then yields 6U exponentially increasing 
(corresponding to U slowing down) with a scale length 

1/2 


L 'v. 


p(0) + 

3 k l 


4 r 2 
P 1 C F 


1 J 


R 


(M 


1) 


1/2 


(3.11) 


where = U./C„ . About the downstream flow, U < C , and (5.10) 

F 1 1 h l 2 F 2 
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yields a damped trailing ICR wave train, as predicted by MacMahon [1968] . 

and Fredricks and Kennel [1968]. 

The above wave train analysis is valid as long as the small ICR 

expansion of the Chew-Goldberger-bow fluid equations is preserved. For 
2 

weak shocks, M - 1 << 1 , (3.11) predicts L » R ; hence the ICR 
F 1 ‘ 
dispersive wave train should resolve the shock structure. If 

M? - 1 > ^- (P^^/p^Cp ) < 1 , L < R + and the above fluid analysis 

breaks down. 

3.2 Discus sion 

A qualitative understanding of the wave train structure for stronger 
shocks is obtainable by considering the linear dispersion relation for 
the perpendicular magnetosonic wave. A collisionless kinetic Vlasov 
treatment of the full electromagnetic perpendicular dispersion relation 
has been performed by Fredricks [1968] ; the result is summarized in Fig. 4. 
The kinetic theory predicts electromagnetic cut-offs (k -*■ 0) and electro- 
static resonances (k ») at harmonics of the ion cyclotron frequency; 
also note the approximate resonance at the lower hybrid frequency. The 
envelopes along the maximum group velocity points of each dispersion curve 
are lines given by w/k = Cp and w/k = Cg . 

Now consider, as in Section 2.0, the predictions of the fast shock 
wave train structure implied by Fig. 4. First note that no upstream wave 
trains are possible since, at the intersections of with the dispersion 

curves, 3w/9k = 0 , and energy would be converted into, not out of, the 
shock. . The downstream flow speed intersects the curves at all cyclotron 
harmonics, and here the group velocity condition for trailing wave trains 
is satisfied. lienee, naively, the kinetic theory predicts a wave train 



- 18 - 


v/ith a multiplicity of discrete oscillation lengths corresponding to each 
harmonic of ft . 

-r 

Several difficulties with the above arguments are apparent. First, 
a finite Coulomb or anomalous collision frequency, which mast be present- 
in the downstream shock flow, would tend to destroy to fine scale ion 
harmonic structure; hence collisionless kinetic theory probably poorly 
approximates the fluid behavior of the plasma. Second, consider a pulse 
with thickness of order R + or less which is undergoing steepening along 
the approximately non-dispersive part, co/k = Cp , of the dispersion curves 
in Fig. 4. Since the ions traverse the pulse width only once (U > C ) , 
the nonlinearly excited waves associated with steepening do not sense the 
full ion gyro-coupling to the magnetic field. To these waves the ion 
orbits appear almost as straight lines; hence ion harmonic dispersive 
structure given by the kinetic theory for infinite plane waves cannot 
effectively limit the steepening of narrow compression pulses. The waves 
genei'.ated by steepening should approximately obey the two fluid dispersion 
relation of Stringer [1963] and Formisano and Kennel [1969] , and remain 
non-dispersive until electron inertia slows down the fast wave at frequencies 
near the lower hybrid. From Section 2.0, the shock structure is a trailing 
C/ojp wave train. Note that the maximum group velocity envelopes including 
the dispersion of electron inertia at the lower hybrid frequency in Fig. 4 . 
approximates the two fluid dispersion curves in Figs. 2a and 2b. 

In summary, for weak perpendicular fast shocks Mp -1 « ~ (Pj^ + /pCp^) , 
L » R + , and the ICR trailing wave train derived from the Chew-Goldberger- 
Low hydromagnetic equations with first order ICR corrections probably 
describes the shock structure reasonably well. For stronger shocks, ICR 
dispersion fails to limit steepening; steepening continues until a C/o)p 
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or possibly X wave train is formed. Finally, even though only 
perpendicular propagation has been considered, the ineffectiveness or 
weakness of ICR dispersion to limit strong fast shock steepening 
probably also holds for oblique strong shocks. Hence in the following 
sections ICR dispersion will be neglected. 
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4.0 Obliqu e Fast Shocks 

If the restriction of perpendicular propagation is relaxed, ion 
inertia dispersion becomes competitive with electron inertia for angles 
0 > tt/ 2 - Al /M . The Mach number dependence of this transition is 
discussed in Section 4.1 using the linear dispersion relation. Section 
4.2 reconsiders the ion acoustic subshock for oblique propagation; 
high-6 and kinetic theory effects are also discussed. 


4.1 Oblique Mag netic Fast Shocks 

Since the wave train differential equation for the whistler shock 
has been extensively analyzed by several authors [see Sagdeev, 1966; 
Cavalierc and Englemann, 1967], the two fluid linear whistler dispersion 
relation will be employed here to determine the' wave train structure. 

The strong shock limit is stressed. 

If | / 00 p « 1 , the phase speed of the fast ’wave is approximately 

given by [Stringer, 1963; Formisano and Kennel, 1969], 


w 


k 2 C 2 


w. 


2 2 
C A cos 0 
A 


,, , 2.,2 . 2.2 , . 2 J 2 . 2 

(1 + k C./Wp ) 1 + k C /u) p 


r 2 . 2 fi 
C sm 0 

D 

7 ~ 7 

1 ♦ k 2 *J 


+ c. 


T+ 

1 + 


1 


2 2 

(1 + k Z X* ) 


(4.1) 


where 0 is the angle between the flow velocity and the magnetic field 
direction. Figure 5 is a sketch of w/k against k ; intersections of 
o/k with possible upstream, and , and downstream, and , 
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flow speeds are included. In order to clearly separate magnetic and 

+ 

temperature effects in (4.1), 3 is taken somewhat less than unity so 

that in the range k = 0 to kX^ < 1 , the first two terms in (4.1) 
dominate the dispersion relation. Note that since for strong shocks ICR 
dispersion was neglected in (4.1), 3 + is not restricted to be 

t. 

< M /M + as in the calculations of Cavaliere and Englemann [1967] . 

If 0 > tt/2 - M /M + , ion inertia dispersion, C/Wp , in the first 


term of (4.1) exceeds that of electron inertia and increases the fast 
wave phase speed above the hydromagnetic speed, Cp . At the first 
intersection of with the dispersion curve, the whistler stands in 

the upstream flow and forms a leading wave train. Note from (4.1) that 
3w/3k > 'w/ k if kC/o)p « 1 ; hence wave energy propagates upstream 
out of the shock. The oscillation length, found by setting w/k = , 

neglecting C/iOp terms, and solving (4.1) for k , is given by 

COS0, 


2v 


(M? -1) 1/2 W P 

+ 


(4.2) 


where M = U /C , the Alfven Mach number. The magnetic wave train 
A 1 1 A 1 

is also sketched in Fig. 5. 

From (5.1) the maximum whistler phase velocity is w/k = (C cos0)/2 , 

3. 

2 1/2 

attained when kC/w = 1 ; C = (B /4 itNM ) . Therefore the maximum 

r 3. — 

upstream Alfven Mach number for which the whistler forms a leading wave 

train is M = Af~7M~ (COS0./2) ; the scale length is L 'u 2 tt (C/w p ) . 

A j ^ -*■ * 

For Uj > (C a cos6p /2 , the only wave train possible is at the 

interaction of the downstream flow speed. (The second Uj intersection 

is disallowed since, 3w/3k < w/k , energy would be- blown back into the 

shock.) If IL > C c , which probably restricts M < 2 to 3 or 
2 S 2 A 1 
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cos 6^ < (2 to 3) x 2M7iT , the linear dispersion relation predicts a 

trailing magnetic wave train. The scale length, obtained by solving (4.1) 

for k and taking the root of the bi-quadratic for U < (C cos6 7 )/2 , 

a 2 


is approximately 



(4.3) 


or L ~ 7i (C/ uip ) . 

For stronger shocks, IJ. > (C cosG 1 )/2 > (2 to 3) x C, , U 0 < C„ . 

^ 1 3-j i A - £ b 2 

Furthermore, even if U, < (C cos6,)/2 so that a C/ou wave train 
’ 1 a, 1 P 

1 + 

leads, < Cg can occur. In previous magnetic wave train analysis 
crossing of the downstream sonic point, corresponds to what has been called 
the breaking of the wave train [Cavaliere and Hnglemann, 1967; see also 
Section 5.0]. Note that the sonic point can occur locally in the wave 
train and need not always occur at the asymptotic downstream flow; nonlinear 
overshoots can also give U = Cg locally. A collisional or turbulent ion 
viscosity, or the ion acoustic subshock wave train is required for a 
complete shock transition. However, since reduction of U to Cg must 
be accomplished by magnetic and density compression, the magnetic shock 
structure should be approximately that predicted by wave trains. For the 

C/Wp wave train, many leading oscillations are probably necessary for 

+ 

U < C c to occur locally in the shock (see Fig. 5). In the trailing wave 
train U probably decreases through Cg in the leading edge (see Fig. 2). 
The initial magnetic gradient should possess a C/u)p characteristic 
length [Sagdeev, 1966]. Since the ions are decoupled from the magnetic 
field at the sonic point, whether or not a complete or even a partial C/wp 
wave train exits downstream is uiclear; further analytical or numerical 
work is necessary to resolve this difficulty. 
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In summary, for low Mach numbers 


M a < v^T/M 
A, + 


cos8^/2 , the shock 


structure is a leading ion inertia wave train. For stronger shocks, a 
dispersion change to a shorter scale length occurs; the wave train now 
trails and has an electron inertia length. If U = Cg anywhere in the 
shock flow, a dissipation (viscous subshock) or dispersion (ion acoustic 
subshock) discontinuity results; structure of magnetic wave trains for 
U < Cg is unclear. 


4.2 Obliq ue Ion Acoustic S ub shock 

Whether or not an ion acoustic subshock occurs for oblique shock 

flows depends on the downstream propagation angle 6^ and the plasma 3 . 

The fast shock evolutionary conditions require that IL C. - C cos0„ , 

l A 2 l 

the intermediate speed [Kan trow its and Petschek, 1966]; hence for B 2 < 1 , 

U 9 < C results only if 0 > cos * (C Q /C ) . For moderately strong 

2 . 2 A 2 
shocks, M„ > 2 to 3 , however, probably exceeds unity so that 

1 

U, < Cg occurs at all propagation angles. 

An additional constraint on the appearance of an ion acoustic subshock 

2 

is that 3^ oi’ 3 2 <(M + /M ) (cos 0/2) , since at higher pressures the 
maximum whistler phase speed w/k = C cos 0/2 < C '; if 0 = tt/2 , this 
constraint becomes B < v4l + /M . Here the dispersion relation (4.1) no 
longer describes the fast wave but becomes the intermediate wave dispersion 
relation if the last two terms are dropped; for 3 > (M + /M ) (cos 0/2) , 
the fast mode is approximately an isotropic ion acoustic wave and remains 
non-dispersive up to ^ 1 [Formisano and Kennel, 1969]. The ion 

acoustic wave train is now just the electrostatic shock derived by Moiseev 
and Sagdeev [1963] in the 3 00 limit. 

A further consideration on the occurrence of the ion acoustic subshock 
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is the relation of U„ to the ion thermal speed. C . Figure 6 for 

w p > |o j cos 8 and Fig. 7 for {Q_jcos0-> are sketches of the 

dispersion curves for the three quasi -hydroma gnetic waves which propagate 
below the electron plasma frequency [Stringer, 1963]; the curves are 
drawn for cos0 < Cg/C^ , < 1 , and C + > Cj . If 0 > 1 , Figs. 6 

t- 

and 7 are modified by moving the sonic line C., closer to the C,. line, 
thus reducing the region of C/oJp dominance in the dispersion relation. 

Note that in the vicinity of kA^ 1 , if a>p > |cos0 , the fast 

+ 

wave first speeds up to Cg and then slows to C + . If |L! |cos9 > to,., , 

the intermediate wave, which is an isotropic sound wave for kC/tOp » 1 

+ 

and kAp << 1 , slows to C + and then intersects the fast mode near 
w r v |cos0 ; the fast wave passes non-dispersively through Cg and then 
speeds up to C + 

Two points merit discussion. First, for the fast shock to steepen 
to an ion acoustic subshock, the linear fast wave dispersion relation 
must possess a non-dispersive region with co/k = Cg . Therefore from 

1/2 ' 

Fig. 6, only shocks in plasmas with uu > [ft |cos0 or > |Q + ft J 

* + - + 

if 0 = ir/2 have ion acoustic subshocks. Second, the downstream flow 
speed for Mp » 1 is constrained by the Rankine-Hugoniot relations to 
be at least V!p , but the temperature jump across the shock has -*• 00 
for M .-*• 00 [Anderson, 1963]; hence U 0 'v C is possible for strong 
shocks. The fluid dispersion curves in Figs. 2, 6, and 7, are non-dispersive 
for w/k = C + , and further shock steepening might be thought possible. 
However, the kinetic theory predicts heavy ion Landau damping when w/k = C + 
[Stix, 1962], so that whether the waves with to/k C + generated by 
nonlinear steepening can propagate is extremely doubtful. Furthermore 
the ion acoustic wave is also ion Landau damped unless T /T >> 1 . 

Therefore the ion acoustic subshock is likely to occur only if 
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co > jQ |cos6 or , C„ - IL » C and T /T + >> I . Strong 

+ 

shock 'flows violating these criteria probably have fully turbulent 
structures, although possibly in conjunction with a magnetic wave train. 

4.3 Discuss ion 

The fast shock wave train structure exhibits 'three dispersion changes 
or discontinuities with increasing shock strength. For 6 = tt/2 and 
M 2 -1 > t (P, (0)+ /PC“) ICR dispersion ceases to inhibit shock steepening, 
and short scale length C/oOp or wave trains develop; if 0 / tt/ 2 

and $ + > 1 so that R + > C/Wp , a magnetic ion inertia wave train is 

+ 

also possible. The whistler magnetic wave train changes from leading with 
C/a)p to a trailing with C/a)p lengths when > /M + /M cos0/2 . 

Finally for Mp > 2 to 3 , U 2 < , and the ion acoustic subshock 

proposed by the Texas group is possible if C*. > U ? » C + , Wp > |fi_|cos9 

|Q + fi_| 1/2 and T"/T + » 1 . 


or 
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5.0 Strong Turb ule nce Modification of Wa ve Train Di spersi on 

In plasmas where the Coulomb collision mean free path greatly exceeds 

even the longest dispersion scale lengths, collisionless plasma turbulence 

must provide the necessary shock dissipation. Collisionless turbulence 

often requires strong gradients, either in Cartesian or velocity space, 

as sources of unstable free energy [Sagdeev and Gal'eev, 1969]. In addition 

to whatever background turbulence exists in the flow, wave train gradients, 

especially in strong shocks, should generate intense turbulent dissipation. 

Wave dispersion properties, and hence wave train scale lengths, however, 

may be modified by strong turbulence. In this section this dissipation 

modification is illustrated, albeit not self-consistently , by a model 

shock flow in which an effective or turbulent electron-ion collision 

1/2 

frequency, > exceeds the lower hybrid frequency \ . 

Note that the usual two fluid equations require v e gf <<: [see Cavaliere 

and Englemann, 1967]. Since ions are now effectively decoupled from the 
magnetic field by collisions, ion inertia and ICR dispersion will not 
limit shock steepening or form wave trains ; hence the shock wave train 
structure should be controlled by electron dispersion. 

To obtain further insight on the collisional suppression of ion 
dispersion, consider the cold plasma oblique fast shock wave train analysis 
of Cavaliere and Englemann [1967]. The ratio of the growth or dissipative 
length, L , of the leading ion inertia wave train to the dispersion 
scale length L was found to be 
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L g (U/V eff ) [cos 2 6/M a 2 (H a 2 -1)3 '(M + /M_) 
L fcose/(M A 2 -l) 1/2 j (C/Up ) 


fi. .. cosG /M 7 h 
lH + - 


V eff M a CM a 2 -1) 1/2 


(5 . 1 ) 


From the discussion of Section 4.1, the maximum value of (cosB v^T + /M ) / 
for ion inertia dispersion to dominate the wave train is 2 \ hence the 

growth length becomes shorter than the dispersion length if v e ff/ftyj ~ 

7 1/7 7 

2/ ( M a -1) • Therefore for sufficiently strong shocks, -1 >> 1 , 

and v e ££ ~ ft j ion inertia scale length effects are suppressed by 

collisions, and the wave train structure must be controlled by electron 

dispersion. Weak shocks under the above conditions are probably adequately 

described by hydromagnetics. ICR dispersion is eliminated if v c ££ > ft + 

since, by the usual arguments leading to the MUD equations, the ion 

pressure tensor will be approximately isotropic. 

The equations describing the strong shock electron dispersion 
structure are the exact conservation relations for the fluid mass, 
momentum, and energy, which, if viscosity is neglected, do not involve 
the dissipation directly, and Ohm's law, which describes the interaction 
of the fluid and magnetic field including dissipation. In two fluid 
theory Ohm's law correct to order M-/M can be written 


K 3jf Cf-f) 


V*B 
— + ~~C~ 


dV 1 «■ 

' M + dt " I 


M - V eff- 

Ne 


(5.2) 


is the fluid velocity, V~ the specific particle velocity, 

J = Ne(V + -V ) the current density, and II the electric field. 

To formulate Ohm's law consistent with v ^ > ftjjj , Eq. (S.2) must 
be ordered. Since electrons, because of their small inertia, carry the 
predominant part of the current, J_ ^ -Ne _V . If the current carrying 
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electron velocity V is of- order ✓M /M_ V + and E ^ V + B/C , the 
terms in (5.2) can be ordered as 


1 . . 1 . 1 V eff 

ft, ,.T “ ‘ - ft f ' ft t • ft, . , 

LH + + LH 


(5.3) 


^rTT v eff 


where T is the time required to go a scale length t L/U . On 
electron inertia scale lengths . t 'w C/o) p U 'w M /ft • ; hence ft r 't AT~/M~ « 1 . 
Therefore from (5,4) and (5.2) consistency requires that over short 
electron dispersion lengths 



(5.5) 

(5.6) • 


Since the mass velocity V_ ^ V , (5.5) predicts that tlie ion motion 
is decoupled from the magnetic field, i.e., the ions obey hydrodynamics. 
Furthermore (5.5) coupled with the momentum equation implies that 

J_x£ 

_ _ Vp‘ = 0 (5.7) 

or that on C/Wp scale lengths, current forces are balanced primarily 

by electron pressure gradients. In (5.6), which replaces the full Ohm's 

law, the only scale length which appears is electron inertia. The 

ordering condition ft^L/U < 1 , however, requires that the flow scale 

- 1/2 

length, L ^ U/ft ^ C/to ^ 0 R ; hence if 8 ^ 1 , electron 

cyclotron radius (ECR) corrections must be included in the wave train 
analysis-.^ 
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The wave train will again be investigated by considering only the 
linear response about the asymptotic stationary flow. To further simplify 
the calculation, the ion pressure tensor is taken isotropic, and the 
electron parallel and perpendicular pressures are assumed equal, both 
probably reasonable due to the high collision frequency; ECR corrections 
to P are retained, however. In addition, only very oblique fast shocks 
are treated so that cos 9 = B /B « 1 , although 6 is not restricted in 

■ TT 

the range - v^1_/M + < 8 < ; consequently U « Cj = C, cosG also holds. 

With these approximations the perturbed fluid equations neglecting 
viscosity and zero order heat flow are 


m B 6 B 

pU <5U + 6P l J + - = 0 

H xx 4 t r 


(5.8) 


ri . B 6 B 

pU oV + <5P (1) - ^ = 0 

y xy 4 ir 


(5.9) 


m B 68 

pU 6V + 6P 1 J = o 

z xz 4 tt 


(5.10) 


pU[ll6U + V 6V ] + y (liSP*- 1 - 1 + P^ 1 - 1 6U) + U5P*- 1 ) + P*- 1 ' 5U 

Z Z c> XX XX 


U B 6B 

, n (l) 1 Z 1 Z 

+ V oP^ ' + -r— 

z xz 4w 


(6q 


,1(1) 






«(1) 


(5.11) 


C 2 i2 ^z C 2 

U ± + — x V 


doB 


O o- v r . = USB + B 6U - B 6V 

_ ,2 2 eff dx z z x z 

0 ) p dx w p 

„2 d 2 6B 


2 d6B 

, - w + ~r~ v C* --T 1 - = U6B. - B fiv 

2 ,2 2 eff dx y x y 

w p dx Up 


U 


(5.12) 


(5.13) 


U l\ = U B z - [8 x 2 (q-B Zj )]/ 4 T,pU 


(5.14) 



u ( 0 )- B d6v 

P ' 7 M 


spd) . 6 pO) * £1- .jl __ 

XX fffl f B d: 


(5.15) 


6P^ 1 - ) 

xy 


n (0) - B ... - 
P v J z d6iJ 

2flTT B~ ~W 


(5.16) 


6P 


(1) _ .ptOJ- B x B 2 dS V> 

XZ 2 Jfi j 


B' 


dx 


2 

, x Ml) x »»UX ~ 5 (P M ') 2 B z d 2 dU' 

(%i + % > * - T ^-72 T “2" 


(5.17) 


(5.18) 


6V 


r dSB 
C z 


y 4rr Ne dx 


(5.19) 


ddB 


6V 


4irNe dx 


(5.20) 


U" = U ; 


6U = 6U 


(5.21) 


In the above equations V = B = 0 at the Rankine-Hugoniot stationary 
points was used. 

Equations (5.8) through (S.21) are straight forwardly reduced to a 
pair of wave train differential equations for <$B z and dB^. except 
for the elimination of 6u from (5.18). Since (5.18) is an ECR term, 
presumed small, and by the assumption of U » Cj , a sufficient 
approximation, in (5.18) only, is to take 6u ~ - (U/B^) oB z . The 
results are 


C 2 

2 

“P 


d 2 6B \ (P (0) " 

£ + £ 

2 2 2 
dx U - C„ 


/P) (B 2 /B 2 ) + (B 2 /4tt P ) d 2 6B 2 

z - R z _£ + 


dx 


w r 


V r e d<SB 
err z 

U dx 


2 2 2 2 
(U -Cp J (U -Cg L ) 

5 9 T‘ 

IT (U-C^) 


6B 


(5.22) 



r 2 d 2 oB ’ r 2 

C y ^ C 

~~ 2 ~ ~ 2 " + T 

0J p dx (Op 


v __ ~d6B 
eff y 

U 


2 2 
U - Cj 2 


dx 


6B 


y 


(5.23) 


2 fOl 2 

where R_ = P l ' /2|ft | NM is the electron cyclotron radius; Cp and 
CgL are the fast and slow hydromagnetic speeds defined by 


'F 

,2 

'SL 


_ 2 ^ 2 

- 

r „ 2 n 2' 

2 1 

- A + C s + 


C . + C c 

A S 

„ 2 „ 2 2„ 
- C. C„ cos 0 
A S 

2 

__ 

2 

k » 


1 1/2 


(5.24) 


Note that 5B and 6B are decoupled, a consequence of the very 
z y 

oblique approximation. 

Solutions of (5.22) describe an exponential rise in at the 

leading edge of the shock and a trailing C/tOp - R wave train, each 
with a scale length given approximately by 


L % 


I (P (0) '/P)(B, 2 /B 2 ) + (B, 2 /4TTP) 2 


1 1/2 




u 2 - C s 2 


— R 


[ (C F 2 1 Mp 2 - 1 1 ) / (U 2 -C S ^)] T/2 


(5.25) 


where L is to be evaluated at the upstream or downstream flow conditions, 
respectively. The damping length for the trailing wave train is 


T (P C0) "/P)(B, 2 /B 2 ) + (B, 2 /4np) 2 


“ 1 1/2 . 


2 

“p 


2 2 
u - c s 


R 


L 'v 
D 


(C 2 /w 2 ) (v eff /U) 


(5.26) 


evaluated for downstream conditions. Note that the familiar wave train 
"breaking" occurs when U < since the pseudo-potential term on the 
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right hand side of (5.22) changes sign and L becomes imaginary. Either a 
turbulent viscosity or a Debye length ion acoustic subshock is required to 
further describe the downstream shock structure. Equation (5.23) does not 
yield a complete shock transition for B , and probably just represents the 
damped rotational intermediate wave which is decoupled from the fast wave when 
ion inertia dispersion is neglected. Hence the magnetic shock should be 
primarily compressive. 

The above analysis is undoubtedly a poor approximation to the ion 
dynamics since changes of the ion fluid on either long, C/uip > R + , or 
very short, X^ , scale lengths are not included. The shock structure 
observed in the magnetic field and electron fluid, however, should be 
reasonably well described by the C/tOp - R trailing wave train up to, 
and maybe also beyond, the sonic point. 

In summary, if strong plasma turbulence creates an effective collision 


frequency v e gf > Hyj » only electrons couple to the magnetic field and 
provide magnetic dispersion. Strong fast shocks have a C/u)p - R_ 
trailing magnetic wave train. Although the above calculation was 
restricted to very oblique fast shocks, the qualitative features of the 
wave train structure are probably valid for a larger range of propagation 
angles. Furthermore note that M > /'M/M cosO/2 was not required to 
obtain a C/u^ - scale length wave train. 
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6.0 Discussion 

6 . 1 Summary 

By use of essentially linear arguments the- following estimates on 
the wave train structure of :fast shocks have been derived: ' 

a) Weak perpendicular finite-3 shocks, Mp^-1 < ^ + /pCp^) < 1 

have a trailing wave train. For stronger shocks, probably including 

oblique propagation, ICR dispersion does not limit steepening; wave 
trains for both 3 < 1 and 3 > 1 will not involve R scale lengths. 

b) Oblique magnetic whistler shocks have leading ion inertia wave 

trains if < /M^/M cos6^/2 and trailing electron inertia wave 

trains if M > AP/M cos 9 /2 . 

A 1 

c) When M > 2 to 3 , the Debye length trailing ion acoustic 

1 ' 2 
subshock proposed by the Texas group occurs if 3 < (M /M ) (cos 0/2) and 

io > |ft | cos 9 , or for 0 = "n/2 , 3 < /M /M and u) n > ft... , 

r - + - r Lrl 

- + . * ■ . 

C > U, » C , and T /T >> 1 . As in its dissipation analogue, the 
b 2 2 

viscous subshock, the ion acoustic subshock decouples the fluid from the 
magnetic field and acts primarily to heat the ions without affecting the 
magnetic shock structure. 

d) Strong turbulent dissipation modifies plasma dispersion 

characteristics. For v __ > ft... and strong fast shocks, ions are 

eff LH 

decoupled from the magnetic field and magnetic wave trains trail with 
C/ p - R scale lengths. Here numerical computations [Kellogg, 1964; 
Bardotti, et al., 1966] are necessary to treat the turbulence generated 
by wave train gradients self-consistently . 



The fluid arguments suggesting that strong shocks might have short 
scale length wave trains are probably somewhat oversimplified. In uniform 
plasmas short wavelength oblique linear waves are often heavily Landau 
or cyclotron damped. Therefore in order for a shock to steepen into a 
short dispersion length wave train, the nonlinear excitation or growth 
rate of these modes must exceed the linear, and possibly nonlinear, damping 
rates. Failure to do so may constitute a criterion for the development 
of a fully turbulent strong shock. It should be noted that Debye length 
wave trains for an electrostatic ion acoustic shock have been observed 
in the laboratory [Taylor, et al., 1970]. Furthermore short C/um 
magnetic scale lengths have been observed in the earth's bow shock. 

6.2 The Earth's Bow Shock 

Although laminar low Mach number shocks have been studied in the 
laboratory for many years [Paul, et al., 1965; Kurtmullaev, et al. , 1966; 
Robeson, et al., 1968], high Mach number, high-3 collisionless shocks are 
presently accessible only in space plasmas. The bow shock formed by the 
interaction of the super-magnetosonic solar wind with the earth's 
magnetosphere ha.s recently been probed by the high telemetry rate satellite 
0G0-5 [Fredricks , et al., 1968]. Although highly variable, the shocks 
examined by Fredricks, et al., [1970] are often characterized by: (see 

Fig. 8) 

1. A magnetic field precursor with scale lengths 'v C/u)p ; 

2. a large magnetic jump followed by several oscillations both 
with gradients ^ C./Wp ; 

3. .large amplitude electrostatic turbulence with frequencies r * 

+ 

which maximize at the C/w p gradients in the magnetic field and are 
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thought to be generated by either the ion acoustic or Bunemann current 
instabilities; 

4. rapid proton thermalization occurring in regions where the 
electrostatic wave amplitudes maximize. 

Theoretically, the most puzzling aspect of these shock observations 
is what appears to be a C/ojp trailing wave train in the magnetic field. 

► 4 " 

For typical solar wind flow conditions, M n, 4-10 , 3“ % 1 , 6 f ir/2 , 

the criteria for a C/top dispersion discontinuity discussed in Sections 
4.1 (M > Ai + /M cos 6/ 2) and 3.1 (0 = tt/2 , M ~ 1.5) are generally 

not satisfied, although exceptions almost certainly occur. Hence on the 
basis of phase standing arguments ion inertia should dominate the longest 
s cal e 1 en gth wave train. 

A possible explanation of the C/tOp scale length is the turbulent 
dissipation modification of plasma dispersion characteristics. For the 
current driven electrostatic turbulence observed in the shock precursor 
and wave train structure, Sagdeev [1965] has estimated the effective 
collision frequency as v ^ = o>p /10 . In the bow shock B 'v 10 ^ gauss 

-3 0 

and N n* 10 cm , Wp /lOSTi^ ^ 10 ; hence, by the arguments of Section 5.0, 

+ 

the ions should be decoupled from the magnetic field, and the magnetic 
wave train should be characterized by the electron dispersion lengths 
C/Wp - R . Ion turbulent decoupling might commence in the precursor, 
thus permitting a C/ 00 p - R wave train in the center of the shock, or 

it might be a feature of the solar wind. Recently Forsland [1970] has 
investigated electron heat flow "current" instabilities in the solar wind 
as a means of turbulently heating ions. If the dissipation rate from 

these instabilities is also of order Wp /10., the ion fluid behaves hydro- 

+ 

dynamically in the solar wind; the electron fluid then carries the magnetic 
field and provides the only magnetic wave train dispersion lengths. 
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6.3 Slow .Shocks in Neutral Sheet Flows 

A popular hypothesis in space and cosmological physics is that 
magnetic energy can be rapidly converted into particle and flow energy at 
x-type neutral points [Levy et al . , 1964; Axford et al., 1965; Sturrock , 
1967]. Petschek [1964] and Petschek and Thorne [1967] have proposed a 
hydromagnetic flow in which most of the magnetic field annihilation occurs 
across standing slow shocks bounding the neutral sheet. It is somewhat 
disturbing that although the neutral sheet in the tail of the magnetosphere 
has been probed by many satellites, no direct evidence of the predicted 
slow shocks has been found. Furthermore typical neutral sheet thicknesses 
are quite thin, ^ 600 km [Ness, 1965] or on the order of C/w n or R . 

r + 

A- 

Therefore if slow shocks indeed exist in the tail of the magnetosphere , 
an ordinary two-fluid hydromagnetic model for slow shocks, which would 
predict shock thicknesses C/oOj, or R + , appears to be inadequate to 

■f 

account for the observed thin neutral sheets. 

To obtain rapid magnetic field annihilation in collisionless. slow 

shocks requires a large turbulent resistivity which probably can only be 

provided by ion acoustic or Bunemann current instabilities driven by steep 

magnetic field gradients. A reasonable speculation is that if the 

turbulent collision frequency satisfies v e ff > ^ the slow shock 

structure is dominated by short scale length electron dispersion with the 

ion fluid decoupled from the magnetic field [Section 5.0). If the BCR 

term is- neglected, Eq. (5.22) also describes a low-3 oblique slow shock 

wave train. At the leading edge B decreases exponentially with a 

z 

scale length 


L a, 



M 

__SL 

.2 ; ,i /2 


‘' ! sl _1 I 


C 



( 6 . 1 ) 



where = U/C^ ; downstream the wave train trails with a wavelength 

given by (6.1) evaluated about the downstream flow. Of course, if 

3 ~ 1 , ECR corrections to (6.1) must be included; the ECR term in (5.22) 

2 2 

is incorrect for slow shocks since terms of order B / B have been 

x 

dropped. In the tail of the magnetosphere C/.Wp 'u 5 km so that a 
complete slow shock wave train could easily l'eside within the "thin" 


neutral sheets observed by satellites. 
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Figure Captions 

Figure 1 . - 

Sketch of the shock structure in velocity and magnetic field for a 

collision dominated hydromagnetic perpendicular fast shock. The magnetic 

2 

shock is characterized by the magnetic Reynolds length r m = C /4iTaU , 

where a is the conductivity based on electron-ion collisions. If 

U 9 < C , a viscous subshock occurs in the downstream flow; ion-ion 
b 2 

collisions slow and heat the plasma across the viscous Reynolds length 

r = y/Pb . where y is the coefficient of viscosity, 
c 

Figures 2a and 2b. 

The perpendicular -fast- linear wave dispersion relation, w vs. k , 

is sketched for op > | and for g~ « 1 (Fig. 2a) and 

± + 

3 >1 (Fig. 2b). Possible intersections of the downstream flow speed 

with the dispersion curves are included. Shock steepening occurs in 

2 2 1/2 

locally non-dispersive regions, w/k = (C^ + C g ) = Cp and w/k = C„ 

steepening is limited by dispersively propagating waves with kC/up 1 
•or kX Q 'v* 1 . The magnetic shock is described by a C/up length 
trailing wave train. If < Cg , intersects the dispersion curve 

near kX Q n, 1 ; an electrostatic Debye length ion acoustic wave train 
occurs in the flow velocity and temperature, but probably not in the 
magnetic field. For reasonably strong shocks, Mp^-1 n, 1 or for g^* n, 1 
the appropriate dispersion curve for estimating wave train scale lengths 
is Fig. 2b; note that the region over which C/up dominates wave 
dispersion is reduced. 
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Figures 3a and 5b. 


The ratio of the downstream sound soeed to flow speed, C r /U , as 

b 2 2 

determined by the Rankine-Hugoniot relations for a perpendicular fast 

shock is plotted in 3a against upstream Mach number M with as 

a parameter. In 3b the critical upstream Mach number, M cr j^ > f° r which 

U 2 = Cg is plotted against' . An ion acoustic or viscous subshock 

results when U 0 < C ; . for M > 2.5 , IL < C for all upstream $’s 
2 S 2 F 1 2 S 2 


Figure 4. 

The perpendicular fast wave dispersion relation as determined by 
Fredricks [1968] from kinetic theory using the full Maxwell equations” is 
sketched. Electromagnetic cut-offs (k 0) and electrostatic resonances 
(k -*■ occur at harmonics of the ion cyclotron frequency, n£l + . The 
envelopes of the maximum group velocity points on each harmonic branch 
are given by the lines w/k = C p and w/k = C g . Note the approximate 
resonance at the lower hybrid frequency, , sketched here to lie 

between (n+1) and (n+2)fl + . 

Figure 5. 

A sketch of w/k against k for the oblique whistler dispersion 

relation. Ion inertia dispersion, C/Wp , increases the phase velocity 

+ 

up to the w/k = C cos8/2 ; electron inertia dispersion, C/w , then 

decouples the wave from the magnetic field and decreases w/k to Cg . 
Debye length dispersion, , further reduces w/k to the ion thermal 

speed C + . If the upstream flow speed intersects the dispersion 

curve near kC/w *\» 1 , the magnetic shock structure is a leading C/w p 

+ + 

wave train. For U. > C cos0/2 no waves stand upstream, and the magnetic 

1 E 
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shock has a short scale < C/ai trailing wave train. If U < C_ occurs 
locally somewhere in the shock flow, the velocity and temperature have an 
ion acoustic subshock wave train with X p lengths. 

Figure 6. 

A sketch of the three quasi-hydromagnet ic waves which propagate 

below the electron plasma frequency [see Stringer, 1963]. The dispersion 

curves are drawn for w > j cos8 , Ml, and C c > C T > C 

+ 

Near kX^ % 1 , the fast wave phase speed first increases to Cg and 
then decreases to C + . 

Figure 7. 

Sanie as Fig. 6 except that jft_J cos0 > w p . The intermediate 
wave, w/k = C g for kC/w p » 1 , decreases its phase speed to C + 
near kX p 1 and then intersects the fast branch at to 'v [ft| cosG . 

The fast wave passes non-dispersively through , and then speeds up 
to C . 

r 

Figure 8. 

Data from a single crossing of the earth's lew shock on March 12, 

1968, as observed by OGO-V [Fredricks, et al., 1970], The first graph is 
the total magnetic field j_Bj in units of 10 ^ Gauss. A short wavelength, 
'v C/u)p , magnetic precursor leads the major increase in the field strength; 
a wave train- like series of oscillations with scale lengths again n- C/iOp 
trail in the downstream flow. The second graph is the ion flux, J + , 
from the Lockheed ion spectrometer. Ion thermalization occurs at the 
magnetic jump. The third graph is the AC electric field spectral amplitudes 
plotted against frequency; the numbers correspond to those in the first 
graph of | B [ The maximum electric field amplitudes occur at the 
maximum gradient of j B [ . 




shock with viscous 





downstream 


^2 

i 

— ^e^fi/p u 

O X 


i 



Multiply dispersive perpendicular fast linear 


/5 ± « 


wave 


doj/dY. < w/k x < 5 / A* 


<3cu/c9k < ca/k v o*?/ / 


c/V/ ^ / / \ 

/ iy 7 \ 

^\ / u i \ 

\ / 

\ Ion acoustic Xq 

\ wave train 


H c/w p_ 


c/cop magnetic wave 
train 


kXn ~1 


ngure 2 









A 


Oblique whistler mode dispersion relation 
co/k vs k 




Figure 5 


ft 


/ 

/ 



Figure 6 


\<* 




Figure 7 


23 ll 55 3(r 


23 u 58 1 GO 


Electric Field 
Spectral Amplitudes 



Figure 8 





